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We have studied 4 He films adsorbed in two porous glasses, aerogel and Vycor, using high pre- 
cision torsional oscillator and DC calorimetry techniques. Our investigation focused on the onset 
of superfluidity at low temperatures as the 4 He coverage is increased. Torsional oscillator mea- 
surements of the 4 He-aerogel system were used to determine the superfluid density of films with 
transition temperatures as low as 20 mK. Heat capacity measurements of the 4 He- Vycor system 
probed the excitation spectrum of both non-superfluid and superfluid films for temperatures down 
to 10 mK. Both sets of measurements suggest that the critical coverage for the onset of superfluidity 
corresponds to a mobility edge in the chemical potential, so that the onset transition is the bosonic 
analog of a superconductor-insulator transition. The superfluid density measurements, however, are 
not in agreement with the scaling theory of an onset transition from a gapless, Bose glass phase 
to a superfluid. The heat capacity measurements show that the non-superfluid phase is better 
characterized as an insulator with a gap. 
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I. INTRODUCTION 

For a wide variety of substrates, adsorbed films of 4 He 
thicker than approximately two atomic layers are super- 
fluid at zero temperature. As the temperature is in- 
creased, a transition occurs to a non-superfluid phase at 
a critical temperature T c . Although the details of this, 
transition depend on the topology of the substrate,trcl 
the phase diagrams for films adsorbed on substrates such 
as Mylar, porous glasses, or packed powders are similar. 
In each case, there is a line of transition temperatures in 
the density-temperature plane separating non-superfluid 
from superfiuid films. At T — 0, this line terminates 
at a critical coverage n c , below which superfluidity does 
not occiir. Only cesium substrates, which are not wetted 
by 4 He,H and the atopiically ordered substrates graphite 
and molecular H^p'Q lead to exceptions to this general 
picture. 

Although the existence of a non-zero critical cover- 
age n c has been known since the earliest experiments 
on unsaturated 4 He films, the onset of superfluidity as 
a function of coverage at T = has received much less 
experimental attention than the superfiuid transition at 
T = T c . This has followed in part from the difficulty 
in making thermodynamic measurements on extremely 
thin 4 He films, but also from a conviction that the non- 
superfluid film of coverage n c forms an essentially inert 
pseudo-substrate for the overlying superfiuid film. Most 
aspects of thin-film superfluidity can indeed be under- 
stood without any consideration of the critical coverage 
n c , and the 4 He coverage is often renormalized by taking 
the effective density of the film to be n — n c . For example, 
the superfiuid density of films with transition tempera- 
tures above 200 mK is observed to be proportional to 
n — n c for many substrates.u 

In contrast to the inert-layer model used for the in- 
terpretation of most experiments, there has been con- 
siderable theoretical work suggesting that the onset of 
superfluidity at T = in 4 Hc films is analagous to a 
metal-insulator transition in Fermi systems. This ap- 
proach follows, from a suggestion of Hertz, Fleishman, 
and Aadersona and has been pursued in depth by Fisher 
et al.pQ who exploited the existence of a natural order pa- 
rameter in Bose systems to develop a scaling theory for 
the onset of superfluidity in strongly disordered systems. 
The spirit of this approach is that the onset transition 
is driven by the competition between exchange, which 
favors superfluidity, and the combined effects of disorder 
and the repulsive He-He interactions, which favor local- 
ization. The localized phase in the case of a 4 He film 
adsorbed on a strongly disordered substrate is presumed 
to be a "Bose glass," with a gapless excitation spectrum 
and a correspondingly nonzero compressibility. 

This paper presents superfiuid density measurements 
of 4 He films adsorbed in aerogel glass and heat capacity 
measurements of 4 Hc films adsorbed in Vycor for cov- 
erages close to the critical coverage n c and temperatures 



down to 10 mK. Our measurements demonstrate that the 
picture of a superfiuid- i nsul ator transition developed by 
Fisher et al. and otherstiHl3 is qualitatively appropriate. 
The heat capacity measurements, for example, indicate 
that a non-superfluid film is far from inert. In fact, the 
density of states of the non-superfluid film has a charac- 
teristic energy which vanishes as the onset transition is 
approached from below in coverage. We do not, however, 
find good quantitative agreement with the scaling pre- 
dictions of Fisher et al. for the temperature dependence 
of the heat capacity and the coverage dependence of the 
superfiuid density and the superfiuid transition tempera- 
ture. We find that the non-superfluid film is best charac- 
terized as an insulator with a gap as opposed to a Bose 
glass. The scaling behavior of thin superfiuid films is 
complicated by both the large energy scale characteristic 
of physical adsorption and the existence of the thermo- 
dynamic superfiuid transition at non-zero temperatures. 
The thermodynamic transition appears to be the dom- 
inant critical point, at least over the temperature and 
coverage range of our study. 

II. BACKGROUND 

In Fig. |l| we show a generic phase diagram for an unsat- 
urated 4 Hc-film in the temperature-density plane. The 
film is superfiuid at T = if the coverage exceeds a criti- 
cal coverage n c , which is typically 1.5 to 2 atomic layers. 
(All of the substrates discussed in this paper are dis- 
ordered and do not support layer-by-layer growth. An 
"atomic layer" thus corresponds to an average coverage 
of one layer.) A phase boundary, corresponding to the 
superfiuid transition temperature T c {n), extends upward 
from the point (n — n c , T = 0). The film is non-superfluid 
above this line. The phase diagram of Fig. 1 is represen- 
tative of the aerogel and Vycor systems discussed in this 
paper as well as a variety of other disordered substrates. 
This paper is devoted to three questions regarding this 
general phase diagram. First, how can we better under- 
stand the onset transition at n = n c : does it have scaling 
properties like conventional second-order phase transi- 
tions? Second, what is the most appropriate description 
of non-superfluid films? Finally, what is our understand- 
ing of the thinnest superfiuid films, within a few percent 
of the critical coverage n c 7 

A. Bose Insulators 

One view of the onset transition is that the 4 He film 
is solid at densities below n c and thus cannot support 
superfluidity. This viewpoint is equivalent to drawing a 
vertical line in Fig. [l] at n = n c and therefore makes a 
distinction between the "solid" part of the film and the 
"liquid" overlayer, which is superfiuid at low tempera- 
ture. Although this approach, often called the "inert- 
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layer model," is. sufficient for interpreting most superflow 
measurements,!! it eschews particle statistics by ignoring 
exchange between the solid layer and the overlying su- 
pcrfluid film. In its most orthodox form, the inert-layer 
model predicts that the various thermodynamic proper- 
tics of the supcrfluid should be extensive with respect 
to n — ?i c , where n is the total density of the film. In 
practice, this does not appear to be the case. For exam- 
ple, the superfluid density p s is generally not linear in 
n — n c for coverages very close to the onset of superflu- 
idity. This observation has to led to some variations on 
the inert-layer model, including the suggestion that the 
compression of the film as it becomes thiekcr leads to an 
increase in the density of the inert layer .Ej 

A more rigorous approach is difficult, and, although 
outlined almost twenty years agoja has been pursued in 
detail only in the last several years. It is quite natural to 
draw an analogy between the metal-insulator transition 
in Fermi systems and the onset transition in a 4 He film, 
which is a Bose system. In both cases, one can argue that 
the onset of metallic (or superfluid) behavior corresponds 
to the appearance of extended states spanning the system 
above a critical density n c . Unfortunately, the theoreti- 
cal machinery developed for Fermi systems is not readily 
adaptable to the Bose case because of the strong repul- 
sive interactions. There is no meaningful non- interacting 
limit, such as pure Anderson localization, for a disordered 
Bose system. 

It is possible, however, to consider a pure but strongly 
interacting system and then gradually turn on the disor- 
der. The simplest representation of the interactions is a 
hard-core repulsion V that prevents two atoms from oc- 
cupying the same site. The effects of quantum exchange 
are introduced via a hopping matrix element J between 
different sites. This Bose-Hubbard model has been stud- 
ied for one-dimpasional (ID) and two-dimensional (2D) 
lattice systems.E3 In both cases, the system is "super- 
fluid" (meaning a many-body extended state exists) for 
any density as long as the exchange energy J is greater 
than a critical value J c . For weak but non-zero exchange, 
the system will be superfluid as long as there are sites 
free to accomodate hopping. This condition fails only at 
commensurate densities, at which exactly n atoms oc- 
cupy each site. In this case, the system is localized by 
the on-site repulsion and is entirely analogous to a Mott 
insulator in the fermionic case. A gap for particle-hole 
excitations exists, corresponding to the energy needed to 
overcome the on-site repulsion when the next atom is 
added to the system. 

Disorder, which can be thought of in the simplest,case 
as a random fluctuation in the on-site repulsion Vi^i has 
two effects. First, it leads to localization at incommen- 
surate densities since atoms will prefer to occupy sites 
with the lowest on-site repulsion. This can be expected 
to destroy superfluidity in the case of sufficiently weak 
exchange. In addition, however, the disorder destroys 
the Mott insulator at commensurate densities, since the 
distribution in V implies that the penalty for double oc- 



cupancy will be reduced at some sites, thus softening the 
gap. As the strength of disorder is increased further, 
the gap closes even for small exchange, and the Mott in- 
sulator disappears completely. Although the system is 
now gapless, it is still localized (by disorder) for small 
exchange strengths n| T-his gapless Bose insulator is often 
called a Bose glass 

B. Bose Glass to Superfluid Transition 

In the case of porous glasses, the adsorption poten- 
tial is strongly disordered on the energy scale (~ I K) 
characteristic of the exchange in a superfluid 4 He film. 
In this case, the onset of superfluidity within the frame- 
work of the Bose-Hubbard model is most likely to occur 
from the Bose glass phase. A characteristic phase dia- 
gram in the plane of exchange strength J and chemical 
potential fi for the case of strong disorder is shown in 
Fig. 0. In a real experiment, the chemical potential is 
tuned by changing the density, and one follows a nearly 
vertical line in Fig. || as the 4 He coverage is increased. 
The onset of superfluidity occurs at the point \x = fi c . 

The situation sketched in Fig. |^ is reminiscent of a con- 
ventional second order phase transition, with the chemi- 
cal potential fi assuming the role of temperature. Fisher 
et al.cQ have exploited this analogy to construct a scaling 
theory of the Bose glass to superfluid transition, based on 
the ansatz that the correlation length £ for fluctuations 
in the phase of the order parameter diverges at the onset 
of superfluidity as 

t~8- v , (I) 

where the reduced chemical potential 5 = /J, — fi c 
parametrizes the distance to the critical point. Since the 
number density and the phase of the order paramater are 
conjugate variables, phase fluctuations are linked to den- 
sity fluctuations, which propagate at the speed of sound, 
introducing a time scale as well as the length scale £ into 
the problem. Fisher et al. define the characteristic fre- 
quency f2, which is assumed to vanish at the onset of 
superfluidity according to the power law 

n~c z , (2) 

where z is the dynamic exponent. The existence of the 
dynamic exponent reflects the quantum mechanical ori- 
gin of the fluctuations in this problem, as opposed to the 
thermal fluctuations at an ordinary critical point. 

Although rigorous arguments are presented by Fisher 
et at, we note here that conventional hyperscaling will 
hold if we assume that the correlation volume has a time 
dimension (l/fi) as well as a spatial part ( d , where d 
is the ordinary spatial dimensionality. This allows us to 
write down an asymptotic form for the singular part f s 
of the free energy density at the onset transition: 

f s {T = 0)~5( d+z >. (3) 
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An asymptotic form for the superfluid density p s on the 
superfluid side of the transition follows from finite size 
scaling x-i 



p.(T = 0)~^, 



C = (d- 



2)v. 



The singular part of the compressibility, 



-d 2 f 
dp? 



(4) 



(5) 



is analogous to the singular part of the specific heat at a 
conventional phase transition. 

In order to make experimental predictions from these 
scaling laws, one needs to know the dynamic exponent 
z. Fisher et al. prove that z = d for the case of a 
Bose glass to superfluid transition. The proof depends 
critically on the Bose glass having a non-zero compress- 
ibility, or, equivalently, a gapless low-energy excitation 
spectrum. 

The theory up to this point makes predictions for the 
scaling properties of the superfluid density only at T = 0. 
Fisher et al. have extended their theory to thermody- 
namic properties at non-zero temperature using finite size 
scaling in 1/T. In implementing this argument, they as- 
sume that the only energy scale in the problem is set by 
the characteristic frequency Vl. Temperature will thus ap- 
pear in the asymptotic free energy in the form kgT/hD,. 
Taking fcs = Ti = 1, the scaling- form for the free energy 
becomes 



f s (T) ~ <H d +*>/(T/fi). 



(6) 



The superfluid transition (at non-zero temperature) oc- 
curs at a universal value u c of the argument T/Q, so that 



(7) 



In order to reduce any of the scaling laws into forms 
testable by experiment, we need to relate <5 = p — p c 
to the density n, which can be tuned in the laboratory. 
Assuming z — d and using the inequality v > 2/d for. 
the correlation length exponent in disordered systems ,E£I 
it can be shown that 



n— n c — p — fi c 



(8) 



which means that 5 can be replaced by n — n c in all of the 
above scaling laws. After doing so, we obtain the scaling 
laws of Fisher et al.rQ 

p s (0) ~ (n-n c ) f , C = v{d + z-2)- (9) 
T c ~{n-n c ) w , w = zv- (10) 
T e ~\p.(p)]*,x=z/(d + z-2), (11) 

where p s (0) is the superfluid density at T = 0. The 
torsional oscillator experiment discussed below addresses 
the validity of these scaling laws for 4 He films adsorbed 
in aerogel glass. 



1. Heat Capacity 

We now consider the low-temperature heat capacity 
on both sides of the onset transition. Technically, the 
only requirement for the Bose glass is that the excitation 
spectrum be gapless. This admits almost any power-law 
dependence of C(T), although for a true glass, i.e. con- 
stant density of states, C(T) oc T. On the superfluid side 
of the onset transition, the dominant low-energy excita- 
tions are long-wavelength third sound modes, leading to 
a Debye specific heat, C(T) ~ T d for T < T c . 

At onset (n = n c ), the specific heat can be inferred 
from the scaling form for the free energy (Eq. ^) . In order 
for f s (T) to remain non-zero and finite at the transition, 
the temperature dependent part /(T/fi) must assume 
the form 



/(T/n) = f{T5~ zv ) ~ {T8- zu y 



(12) 



since the sum of all the powers of 8 in Eq. ^ must equal 
zero. At onset, 



C(T) = T 



Q T 2 



(5 = 0,T -►()). (13) 



Note that if z = d, as is believed to be the case for a 
Bose glass to superfluid transition, lim,T->QC{T) ~ T at 
onset. The temperature dependence at onset is thus the 
same as in the Bose glass phase. 



2. The Inert-Layer Model 

Fisher et a/0 demonstrate that in the mean-field limit, 
their scaling theory reduces to the inert-layer model, in 
which the superfluid film is taken to coexist with a solid- 
like 4 He underlayer. In this case, exchange between the 
superfluid film and the localized layer is assumed to be 
negligible, and the thermodynamics of the system near 
the onset of superfluidity can be deduced from consid- 
eration of an ideal Bose gas coexisting with a solid 4 He 
pseudo-substrate. 

The superfluid density in this limit should be propor- 
tional to n — n c , where the critical coverage n c is the 
density of the solid layer. Thus, p s (0) ~ (n — n c )S where 
C = 1 for both d = 2 and d = 3. Two more exponents 
can be inferred directly from the density dependence of 
the Bose-Einstein condensation temperature: 



2nTT 



N 



VC(d/2) 



(14) 



where V is the volume (or area in two dimensions) of 
the system. Strictly, T c = for d = 2 because of the 
logarithmic divergence of the zeta function £(!)■ This 
divergence can be removed through the introduction of a 
short-length cutoff which accounts for the finite thickness 
of the film. We thus find T c ~ (n - n c ) w and T c ~ 
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[p s (0)] a: , where w = x = 1 for d = 2, and w — x = 2/3 
for d = 3. 

We now propose low-temperature forms of the heat 
capacity for the inert-layer model. For n > n c , the heat 
capacity of the ideal Bose gas has the limiting behavior 
C ~ as T -> 0. (The presence of superfluidity re- 
quires a Debye contribution C oc T d , but this vanishes 
more rapidly than the Bose gas contribution in the limit 
T — > 0.) Since the inert-layer model makes no assump- 
tions about the non-superfluid film, it admits essentially 
any form of the heat capacity for densities n < n c . 

3. Summary of the Scaling Laws 

In Table |, we summarize the results of the scaling the- 
ory for both the Bose glass and inert-layer models. Four 
possible cases are considered, corresponding to the Bose 
glass or inert-layer models in either 2 or 3 dimensions. 
In calculating the exponents from Eqs. ^- [n] and [l3[ we 
have used the result z =r-d as well as the correlation 
length inequality v > 2/dx3 

III. EXPERIMENTAL DETAILS 

Two experiments will be discussed in this paper. The 
first is a systematic study of the superfluid density of 
4 He films adsorbed in 91% porosity aerogel glass for tem- 
peratures between 10 mK and 200 mK. The second ex- 
periment is a high-resolution heat-capacity study of He 
films adsorbed in Vycor glass for temperatures from 10 
to 200 mK. The heat capacity measurements were con- 
ducted for coverages on both sides of the onset transition. 
Brief accouJits_of these experiments have been published 
previouslyEJlla 

A. Substrates 

Vycor and aerogel are both porous glasses. VycorEl 
is prepared by leaching out one phase of a borosilicate 
glass after spinodal decomposition. The resulting matrix 
has a distribution of pore diameters with a maximum at 
approximately-,70 A and a half-width of approximately 
20 A. AerogelEj is prepared by kinetic aggregation of sil- 
ica particles in a gel. After hypercritical drying to remove 
the solvent, the resulting substrate consists of strands of 
silica approximately 50-100 A in diameter. Unlike Vycor, 
there is no characteristic pore size, but the structure is 
correlated on length scales up to 600 A for the sample 
studied in this work.Q In spite of the structural differ- 
ences, the superfluid densities of thin 4 He films adsorbed 
in Vycor and 91% porosity aerogel near T c are remark- 
ably similar, although the heat capacity singularity at 
the superfluid transition seen in Vycor does not appear 



m aerog elB As discussed in Ref. the correlated disor- 
der of the substrate is apparently irrelevant for very thin 
films of superfluid 4 He adsorbed on either substrate, since 
the superfluid correlation length £ exceeds all structural 
length scales. 

B. Superfluid Density Measurements 

The superfluid density of 4 He films adsorbed in 91% 
porosity aerogeL-was measured using the torsional oscil- 
lator technique.Eil The experimental cell was the same as 
that described in Ref. f|, which discusses the experimental 
aspects of the superfluid density measurements in detail. 
Several important parameters are listed in Table 0. All 
of the torsional oscillator data in this paper will be pre- 
sented in terms of the superfluid period shift AP defined 
below. These data can be converted to a superfluid mass 
by dividing by the superfluid mass sensitivity of the cell, 
0.022 nsec/pg. 

The measurements under discussion here can be di- 
vided into two classes. The first comprises tempera- 
ture sweeps conducted in a manner identical to that de- 
scribed in Ref. [|. For each 4 He coverage, the resonant 
period P and the amplitude of the torsional oscillator 
were measured as a function of temperature. To deter- 
mine the superfluid period shift AP, we first subtracted 
the temperature-dependent background of the empty cell 
from the raw data. The data above the superfluid tran- 
sition temperature T c were fitted to a constant. The pe- 
riod data were then subtracted from this baseline, giving 
the superfluid period shift AP(T). The period shift at 
T = 0, AP(Q), was determined by fitting AP(T) to a 
horizontal line at the lowest temperatures. Uncertainties 
due to the background subtraction and the fact that our 
lowest reliable measurements could be made only above 
10 mK limited the precision of our AP(0) measurements 
to about ±0.1 nsec. 

The superfluid transition temperature T c was deter- 
mined with a typical uncertainty of ±0.2 mK by locat- 
ing the break in AP(T) at the transition. The most 
precise way to do this was to drift slowly (dT/dt ~ 
0.2 — 0.3 mK/hr) through the transition while record- 
ing the resonant period and amplitude of the torsional 
oscillator. The result of such a drift measurement for 
n = 35.7 /imoles/m 2 is shown in Fig. ^|, for which 
T c = 38.6 ± 0.2 mK. These data also show two dips 
in the amplitude. The large feature, near 38 mK, is 
certainly a sound resonance since the period, which is 
measured using a phase-locked loop, shows an s-shaped 
distortion. The smaller peak, which seems to coincide 
with the superfluid transition, may be due to a cascade 
of resonances. We do not know why the resonance asso- 
ciated with the large peak is so prominent, since a simple 
calculation shows that it does not correspond to a low- 
order third sound resonance. For the thinnest films, the 
resonance could be seen even when AP(T) was almost 
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obscured by noise and uncertainties in the background 
subtraction. 

For films with transition temperatures above 50 mK, 
we could also determine T c by fitting the period shift data 
near the transition to a power law AP(T) = At c , where 
t = 1 — T/T c .a The T c determined in this fashion was typ- 
ically 0.2-0.3 mK below the temperature at which AP(T) 
was non-zero within the resolution of our measurements. 
In this sense, the superfluid transition for films on aero- 
gel is considerably sharper than that observed in rods of 
Vycor,ErE3 where a "foot" in p s extends 1-2 mK above 
the superfluid T c as determined from power-law fits.E3 

C. Heat Capacity Measurements 

The calorimeter, shown in Fig. ^, was constructed from 
thin- walled silver. Fourteen disks of thickness 0.3-0.6 mm 
were cut from a slab of Vycor by first coring the slab 
with a slurry cutter and then slicing off the disks using 
a diamond-impregnated wire saw. Cutting from a slab 
with the leaching plane parallel to the surface produces 
a more homogenous sample than is obtained from rods 
of Vycor, in which the boron-rich phase is leached from 
the perimeter .E2l After cutting, the disks were cleaned in 
hydrogen peroxide and then baked in vacuum. A silver 
coating, ~0.5 /im thick, was then evaporated on both 
sides of each disk. The edges were left uncoated to allow 
the 4 He to penetrate the sample. Each disk was then 
glued to two 0.025 mm thick Ag foils, which had previ- 
ously been diffusion-welded to the cap of the calorimeter. 
The calorimeter was sealed with Emerson Cumming Sty- 
cast 1266 epoxy. A Pt-W wire heater was wound around 
the body of the calorimeter, and a 1 kil RuC>2 resistor 
was glued to the cap of the cell. The calorimeter was 
supported by four Vespel posts and was attached to the 
low-temperature stage of a PrNis demagnetization cryo- 
stat. The thermal contact to the stage was controlled 
using a tin-wire superconducting heat switch. 

A low thermal mass DC magnetization thermometer, 
shown in Fig. was essential for measurements below 
20 mK. The actual sensor in this thermometer is a minm^ 
turized version of the design of Greywall and Busch.o 
2.5 mgrams of cerous magnesium nitrate (CMN) pow- 
der, consisting of grains less than 75 microns in diame- 
ter, was mixed with an equal volume of 2-4 micron silver 
powder and then pressed between eleven 0.025 mm thick 
silver foils. Before pressing with the powder, the foils 
were annealed, perforated to enhance the adhesion of the 
Ag-CMN powder, and their ends were diffusion welded 
together. After pressing, the assembly was potted in Sty- 
cast 1266 epoxy and machined into a cylinder 3 mm in 
diameter. A 2 /iH NbTi wire coil, which was part of a su- 
perconducting flux transformer, was then wound around 
the cylinder, and the foils were attached to a Ag wire 
using a crimp joint. This assembly was then repotted in 
epoxy and glued to a Vespel post, which was attached to 



the low temperature stage of the refrigerator. Two Nb 
shields were used to trap a small (25-50 Gauss) magnetic 
field and for magnetic shielding. The wires for the flux 
transformer were twisted and then threaded through a 
Nb capillary to a junction box on the still of the dilu- 
tion refrigerator. The capillary was filled with silicon oil 
in order to prevent any vibrations of the twisted pair. 
The remainder of the flux transformer passed through a 
feedthrough into the bath space and was-attached to the 
input coil of a commercial DC SQUID £3 

The operation of the thermometer is identical in prin- 
ciple to the high-resolution thermometer developed for A- 
point experiments by Lipa and Chui.c!3 Sppe details have 
been published previously by our group .til The magneti- 
zation is measured by monitoring the current in the flux 
transformer, which changes so as to keep the total flux 
enclosed by the superconducting loop constant. The sen- 
sitivity of the thermometer is determined by the applied 
field. The price paid for the sensitivity of the DC mag- 
netization thermometer is the need to count flux, since 
the dynamic range of the flux-locked loop is limited, and 
so the loop must therefore be reset periodically. In prac- 
tice, we regularly exceeded the slew rate of the SQUID 
electronics and the flux count therefore became mean- 
ingless. This inconvenience was ameliorated in our case 
by the availability of the heat switch on the calorimeter, 
which could be closed to check the magnetic thermome- 
ter against a 3 Hc melting curve thermometer on the stage 
at any point during the experiment. 

The heat capacity measurements were conducted us- 
ing the adiabatic calorimetry technique. A current pulse 
was applied to the heater, and the resulting thermal re- 
sponse of the cell was recorded. Two response curves, 
recorded at temperatures of 102.1 and 18.6 mK for a 
coverage of 27.1 /xmoles/m 2 , are shown in Fig. ^. These 
data are representative of two classes of response curves 
that we found. Above 80 mK, the temperature of the 
cell decayed exponentially after the pulse, and the tem- 
perature step AT was determined simply by extrapolat- 
ing the exponential decay back to the midpoint of the 
heat pulse. At lower temperatures, the decay became 
non-exponential, as can be seen in the lower panel of 
Fig. |(| Since we could fit the long-time part of the decay 
to an exponential, we first believed that the short-time 
behavior reflected a thermal overshoot due to poor ther- 
mal contact between the 4 He film and the calorimeter. 
If we extrapolated the long-time behavior back to the 
midpoint of the pulse, however, the inferred heat capac- 
ities were clearly too high. This suggested that either 
heat was leaking out of the calorimeter on short time 
scales, hence leading to an overestimate of the heat ca- 
pacity, or that the anomalous decay was due to some 
internal relaxation mechanism, perhaps reflecting a sig- 
nificant heat capacity in poor thermal contact with the 
rest of the calorimeter. The former possibility was un- 
likely since the superconducting heat switch provides its 
best thermal isolation at the lowest temperatures. We 
suspected that the Vycor itself might harbor a thermal 
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reservoir and modeled the thermal relaxation inside the 
cell assuming the existcn.ce of a reservoir that was weakly 
linked to the 4 Hc film.LL3. This provided a reasonable fit 
to the decays, but we found that the thermal conductivity 
to the weakly linked reservoir depended on both cover- 
age and temperature. A similar coverage-dependent cou- 
pling between the 4 He film and the Vycor itself has been 
observed in ultrasound measurements. t3 We decided to 
proceed on the assumption that the 4 He film equilibrated 
rapidly with the thermometers, in which case the temper- 
ature step AT for determining the specific heat should be 
extracted from the response of the calorimeter at short 
time scales. We tested this hypothesis by dosing approx- 
imately 2 monolayers of 3 He into the calorimeter at the 
end of the run. This film had a heat capacity approxi- 
mately 60 times larger than any of the He films in our 
experiment. The 3 He film equilibrated within 120 sec 
after each heat pulse, justifying our assumption about 
the short time-scale behavior of the calorimeter at low 
temperatures. 

The heat capacity for a coverage of 27.6 /xmoles/m 2 
is shown in Fig. |7] after dividing by the temperature T. 
These data, taken with the magnetization thermometer, 
show a cusp at 20r«iK, which is a signature of the su- 
perfluid transition.c3E3 (Since we did not take a BET 
isotherm for this sample, the surface area was determined 
by comparing our measured superfluid transition temper- 
atures-|as a function of 4 He dosage with those of Crooker 
et al.c3) The ordinary low-temperature limit of our mea- 
surements was 10 mK. This was determined by the in- 
ternal heat leak of the calorimeter, which decayed from 
approximately 0.8 to 0.4 nW over the course of three 
months. This large heat leak, which we suspect was due 
to the Vycor in the calorimeter, required us to run with 
the heat switch partially open at the lowest temperatures. 
Obtaining a proper tuning of the heat switch was very te- 
dious, and this turned out to be the primary impediment 
to using the calorimeter below 10 mK. 



IV. RESULTS 

A. Superfluid Density 

In Fig. H, we show the period shift data for twelve 
coverages of 4 He adsorbed on aerogel. The superfluid 
period shift at T = 0, AP(0), was obtained by fitting the 
low-temperature part of each curve to a horizontal line. 
The transition temperature T c was determined for each 
coverage using the drift technique described above. In 
Fig. ||, we show T c and AP(0) as a function of coverage. 
The curves show fits to power laws 



where we find w = 1.59 ± 0.06 and C = 1-64 ± 0.04. 
The critical coverage n c — 33.89 ± 0.08 /imoles/m 2 was 
determined from the fit of T c versus coverage and then 
held fixed for the fit of the period shift data. The final 
relation of interest is the dependence of T c on AP(0), 
which is shown in Fig. [10[ In this case, we have fitted 
the data to the power law 



T c = A[AP(0)] x , 



(17) 



where we find x = 0.95 ± 0.02. The measured exponents 
are listed in Table [II along with corresponding values 
we obtain from fits of the data of CrookerE3 for the 4 He- 
Vycor system. 



T c {n) = T (n/n c -iy 



(15) 



We will discuss the exponents of Table |I| in more de- 
tail below. Another point of interest is the size of the 
"critical region" for the onset transition. Since the fits of 
the various power laws are poor (by the standards gen- 
erally applied to critical phenomena), the usual practice 
of looking at deviations from the power-law behavior is 
not particularly instructive. Another approach is simply 
to see how the data near onset compare with those fur- 
ther away from the transition. We do this for AP(0) as 
a function of coverage in Fig. [ll[ The highest coverage 
in this figure has a T c of approximately 1 K. A signifi- 
cant departure from linear behavior is seen only for films 
with AP(0) < 1.5 nsec, or T c < 150 mK. The linear 
dependence of AP(0) on n at higher coverages is consis- 
tent with the inert-layer model of the onset transition, so 
the data of Fig. [TT] indicate only a small coverage regime 
in which fluctuations associated with the onset transi- 
tion might be relevant. The evolution of T c with AP(0) 
is even more striking in this regard, since, as shown in 
Fig. O, the nearly linear behavior observed for thin films 
continues for films with T c 's up to 0.7 K. [This linear de- 
pendence is suggestive of the Kosterlitz-Thouless- Nelson 
(KTN) relation for a two-dimensional superfluidjEil for 
which p s (T-)/T c = 8.73 /imoles m" 2 K-l, where p s {T~) 
is the limiting superfluid density as the transition tem- 
perature is approached from below. For reasons of com- 
parison in Fig. O, we have converted the period shift to 
areal superfluid density using the calibration p(s)/AP = 
1.24 /imoles m _2 nsec _1 . The slope of the KTN line, 
which is drawn dashed in the figure, is about twice that 
of a linear fit to the data of Fig. [l2]. Although we expect 
p s (T~) to be reduced from p s {0) by elementary excita- 
tions and vortex pair screening, the discrepancy between 
the two slopes is not resolved in a full fit— j>f the super- 
fluid density data to the finite size theoryE3 for a super- 
fluid film adsorbed on a cylindrical strand. The best-fit 
coverage Jar each data set that is about half the actual 
coverage^ We do not believe, however, that the failure of 
this naive application of the KTN relation says anything 
about the effective dimensionality of the film, which will 
be discussed further below. 1 



and 



AP(0) = D (n/n c -lf 



(16) 
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B. Heat Capacity 



A. Scaling Theory 



The 4 He-Vycor heat capacity data, divided by tem- 
perature, are shown for ten coverages in Fig. Each 
of the data sets in the lower panel of the figure shows 
a sharp cusp, corresponding to a superfluid transition.ESI 
The superfluid transition temperatures determined from 
the locations of the cusps are shown as closed circles in 
Fig. [l4|. The critical coverage n c for the onset of super- 
fluidity is determined by extrapolating the T c data to 
zero temperature. For temperatures T < T c , the heat 
capacity of the superfluid films is roughly quadratic in 
temperature, while it depends nearly linearly on temper- 
ature for T > T c , with a small quadratic correction that 
increases with coverage. In the vicinity of the superfluid 
transition, thepe, is a small peak in C/T due to critical 
fluctuations.c§E3 The size of the peak decreases with de- 
creasing thickness. As can be seen in Fig. |7| this peak 
can barely be resolved for the thinnest superfluid film (T c 
= 20.5 mK). 

At high temperatures, the heat capacity of non- 
superfluid films, shown divided by temperature in the 
upper panel of Fig. |l3|, depends linearly on temperature, 
but C/T drops rapidly at low temperatures. We define a 
crossover temperature Tb separating these two regimes to 
be the point at which C/T is half of its high-temperature 
value. T B is shown as a function of coverage in Fig. [l4|, 
along with the T c 's measured for the superfluid films. 
This plot shows clearly that Tb vanishes as the onset 
transition is approached from below. (Although our defi- 
nition of Tb is somewhat arbitrary, any reasonable choice 
of a characteristic temperature for the crossover would 
show the same coverage dependence.) For T > 2Tb, 
there is not much qualitative difference between the heat 
capacity of a non-superfluid film and that of a super- 
fluid film for T > T c . This is emphasized in the three- 
dimensional representation of our data shown in Fig. [l^, 
in which there is a wedge-shaped plateau in the vicinity 
of the onset coverage n c — 27 /imoles/m 2 . Above 20 mK, 
the heat capacity isotherms, which are cross-sections of 
the 3D surface, show no feature at n c . The peaks in the 
isotherms become progressively narrower in coverage as 
the temperature is reduced, but we cannot reliably ex- 
trapolate our data to T = 0. 



V. DISCUSSION 

The data introduced above indicate that the concept 
of a superfluid-insulator transition in 4 He films is qualita- 
tively appropriate. Both sets of measurements indicate 
the existence of a critical coverage n c at which the su- 
perfluid transition temperature T c vanishes. The heat 
capacity data are particularly instructive in this regard, 
since there is a characteristic energy ksTs for nonsuper- 
fluid films which also vanishes as n c is approached from 
below in coverage. 



1. p s measurements 



As can be seen by comparing the exponents of Table III 
with the predictions of Table |, the experimentally deter- 
mined exponents for the onset transition are not consis- 
tent with the scaling theory for either d — 2 or d — 3. We 
have also seen in Figs. |ll| and [l2]that the departure from 
the behavior seen at higher coverages (T c > 200 mK) 
is actually quite small. Some curvature is apparent in 
AP(0) at the lowest coverages, but the relation between 
T c and AP(0) remains nearly linear for T c 's from 700 mK 
down to the lowest T c in our study (20 mK). These ob- 
servations lead us to question whether our measurements 
are in the asymptotic regime of reduced coverage n — n c 
in which the scaling theory is valid. To address this ques- 
tion, we return to the asymptotic form of the free energy, 
Eq. This scaling form comprises the product of the free 
energy at T = and a function of T/fl = u c T/T c . This 
implies 



Ps(T) 

Ps(0) 



u(u c T/T c ) 



(18) 



in the asymptotic regime, where u is a universal func- 
tion. Thus, upon scaling the period shift AP by AP(0) 
and the temperature by T c , all of the AP(T) data in 
the asymptotic regime should collapse on to a universal 
curve. As can be seen in Fig. |l^, our data do not scale 
in this fashion. In fact, the scaled data in this figure 
drift continuously toward the origin as the coverage is 
reduced. There is no indication that the data are even 
approaching a universal curve. The scaling plot demon- 
strates explicitly that the properties of superfluid films 
at non-zero temperature are not determined by the sin- 
gle energy scale ft. Therefore, we should not expect the 
scaling laws (Eq. |l^ and |lTJ) for T c as a function of n or 
p s {0) to hold in the coverage range of our experiment. 

Although we have shown explicitly that the scaling the- 
ory does not hold at non-zero temperature (i.e. T > 
20 mK), we must still consider the possibility that it 
holds at T = 0. Only one of the scaling laws, p s (0) ~ 
(n — ri c ) i applies in this case, but the experimental re- 
sult £ = 1.64 ± 0.04 is well below the theoretical bounds 
C > 8/3 for d = 3 or £ > 2 for d = 2. Since, as shown 
in Fig. the data approach the inert-layer prediction 
at higher coverages, it is reasonable to ask if our low- 
coverage data fall in a crossover regime in which quantum 
fluctuations are relevant but not yet dominant. This is 
certainly possible, but it leads one to ask why the critical 
region in coverage should be so small. In order to address 
this question, we convert the coverage to the chemical po- 
tential fi using the Van der Waals relation 



/x(n) 



Q 



(19) 



where a = 27 K (layer) 3 for 4 He on glass. Although the 
chemical potential is difficult to measure in this coverage 
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range, third sound measurements, which determinej-it(ri) 
indirectly, are consistent with the form of Eq. [yj.El As- 
suming a monolayer coverage of 13 /xmoles/m 2 , n c w 2.6 
layers for 4 He-aerogel and n rs 3.0 layers for a film with 
a T c of 200 mK. The chemical potential thus shifts from 
about -1.5 K at onset to approximately -1.0 K for films 
with T c 's near 200 mK. Given that the characteristic ex=. 
change energy for the 4 He films is on the order of 1-2 K,E£| 
a chemical potential shift of 0.5 K cannot be considered 
small. 

The rapid change of fi with coverage also leads to a 
technical problem in applying the scaling theory: the as- 
sumption n — n c ~ fi — /j, c is not correct, even in the nar- 
row coverage range of our experiment. Assuming that fi 
follows the form of Eq. [l^, we have converted coverage to 
chemical potential and have fitted our data to the power 
law 

AP(0)~(^p c ) c . (20) 

In effect, this allows us to eliminate the assumptions im- 
plicit in Eq. |^, with the caveat that we have not actually 
measured the chemical potential. We find an exponent 
£ = 2.2 ±0.1 by this approach. The case for the Bose 
glass is thus improved in d = 2, although the inequality 
of Eq. H is still not satisfied for d — 3. We note, how- 
ever, that the theory cannot be applied self-consistently 
in d = 2, since it assumes a diverging correlation length 
which will eventually exceed the strand size of the aero- 
gel. In this case, topological arguments lead us to expect 
that the onset transition must be three-dimensional in 
character, although we cannot exclude the possibility of 
an observable crossover regime between 2D and 3D be- 
havior. 

We thus find that two conditions of the scaling theory 
do not apply in the case of 4 He adsorbed in aerogel in 
the coverage range of our study. First, there is no single 
energy scale that determines the thermodynamic prop- 
erties at finite temperature. Second, the chemical po- 
tential changes rapidly with coverage, which makes the 
asymptotic regime difficult to achieve and also implies 
that the reduced density n — n c is not simply propor- 
tional to n — fi c . We recall at this point the third (and 
most fundamental) assumption of the theory, which is 
that the non-superfluid phase is compressible. (This as- 
sumption is critical to establishing the relation z = d 
for the dynamic exponent.) Since the compressibility is 
given by dn/dfi, Eq. [l9] suggests that although k will 
be non-zero, it will be considerably lower than for an 
idealized 4 He film such as that considered in the Bose- 
Hubbard modelJi9~li3 In effect, the adsorption potential 
adds a third energy scale to the problem in addition to 
the exchange energy and the interparticle repulsion. The 
strong binding to the substrate makes the film much less 
compressible than it would be if it were "free-standing. "L3 
Even if k is technically non-zero, the film may be suffi- 
ciently incompressible that the softening due to quantum 
fluctuations is negligible. 



2. C(T) measurements 

The heat capacity measurements allow us to test the 
predictions of the scaling theory below the onset of super- 
fluidity. As shown in Fig. [l3|, the salient feature of C(T) 
for coverages below onset is a characteristic temperature 
Tb below which C/T drops rapidly to zero and above 
which it approaches a constant. As shown in Fig. [u], Tb 
vanishes as the onset of superfluidity is approached from 
below. At the onset transition, C(T) is approximately 
linear over the entire temperature range, T > 10 mK, 
covered by our experiment. 

The fact that C (T) is linear in T for n = n c is in appafe 
ent agreement with the scaling theory of Fisher et alS3 
(see Table §). Below onset, however, the data show only 
a single region of linear behavior, although the Bose glass 
model predicts two regions: one at low temperature due 
to the Bose glass and a second at higher temperatures 
due to quantum fluctuations. One possible explanation 
for the single linear region is that the Bose glass contri- 
bution is immeasurably small. This would be consistent 
with the observation made above that the strong Van der 
Waals potential makes the 4 Hc film nearly incompress- 
ible. 

Assuming for the moment that the Bose glass contribu- 
tion to the heat capacity is too small to be measured, we 
turn to the linear behavior that we do observe. At onset 
this is predicted by the scaling theory of Fisher et at, in 
which the contribution to C from quantum fluctuations 
varies linearly with temperature. We note, however, that 
the linear behavior is observed over a very wide tempera- 
ture range, extending up to at least 600 mK, while the su- 
pcrfluid density data for both aerogel and Vycor indicate 
that the critical regime in which quantum fluctuations 
are relevant is confined to coverages with T c 's less than 
200 mK. Although there is no reason why the asymptotic 
critical regions observed for two different thermodynamic 
quantities must be the same, it would be rather unusual 
to see critical behavior so readily in C (T) when all other 
measurements suggest that the critical regime is nearly 
inaccessible. 



B. Activation Model for Non-superfluid Films 

As can be seen in Fig. the heat capacity of non- 
superfluid films drops off very rapidly, faster than T 3 , at 
low temperatures. This behavior is suggestive of a gap 
in the excitation spectrum of the film. The simplest den- 
sity of states that includes a gap is shown in Fig. [o. 
The model, considered originally by Tait and ReppytH 
includes two bands of states. The lower band comprises 
localized states due to the heterogeneous adsorption po- 
tential. Because of the strong repulsive He-He interac- 
tions, each localized state can accommodate only a single 
He atom, leading to an effective Fermi statistics in the 
lower band of Fig. The upper band, which is unoc- 







cupied at T — 0, consists of extended states, which are 
assumed to obey Bose statistics at low densities. This 
model leads to a heat capacity of the form 

C = D(A/2k B T + 2)e- A/2kBT (21) 

for k B T < A. The prefactor D = ^NiN u Ak B , where N t 
and N u are the densities of states in the lower and upper 
bands. The results of fits to Eq. |l]for four non-supcrfluid 
films in the regime T < T B are shown in Fig. nM The 
gap A is shown as a function of coverage in the inset. For 
each coverage, we find a gap of order 5T B . 

Our data are not precise enough to distinguish between 
a true gap in the density of states and a "soft gap" , in 
which a low-energy tail exists in the band of extended 
states. It is possible that a small number of extended 
states exist at energies well below A. In principle, it 
would be possible to place a bound on the number of such 
states by determining how much of the classical entropy 
is not frozen out at low temperatures. Unfortunately, this 
calculation is impossible because our data do not reach a 
limit in which C(T) is independent of temperature. (This 
limit could not be reached without macroscopically pop- 
ulating the vapor phase, which is ignored in our model.) 

There is no a priori reason for our having chosen a 
constant density of states in each of the bands shown in 
Fig. [It]. As far as fitting the data of Fig. [l8] is concerned, 
the critical element of the excitation spectrum is the gap. 
Additional information, however, comes from the high- 
temperature data on both sides of the onset transition. 
In both cases, the heat capacity depends linearly on tem- 
perature in this regime. For n < n c , only the extended 
states contribute to the heat capacity at high tempera- 
ture, and a linear temperature dependence would follow 
from a constant density of states. The only difference for 
superfluid films (n > n c ) is that some extended states 
are occupied at T = 0. For T > T c we would expect 
to see the same temperature dependence as observed for 
films below the critical coverage. Based on Fig. [l5], this 
appears to be the case, since the wedge-shaped plateau 
at high-temperatures extends to both sides of the criti- 
cal coverage. We have also fitted the high-temperature 
(T > 2T B or T > T c ) data to the form C(T) = AT + BT 2 
on both sides of the onset transition. As can be seen 
in Fig. [l9], the linear coefficient A increases slowly with 
coverage, with no apparent feature at the onset transi- 
tion, which is marked in the figure by the dashed line. 
The quadratic contribution, however, increases rapidly 
for coverages above onset, although it remains at least 
a factor of two smaller than the linear heat capacity for 
T < 400 mK. 

The onset of superfluidity in 4 He films adsorbed in Vy- 
cor thus corresponds to the closing of a gap in the exci- 
tation spectrum of the film as the coverage is increased. 
Although this is a Bose insulator to superfluid transition, 
the insulating phase is not gapless, in disagreement with 
the assumptions of the Bose glass model. Although this 
conclusion turns to some degree on whether or not we can 



resolve a small but nonzero density of states at low tem- 
peratures, there is also a significant philosophical distinc- 
tion between the activation model we have employed here 
and the Bose glass approach. In the Bose glass model, 
superfluidity is destroyed at n — n c by order parameter 
phase fluctuations. The Bose condensate (essentially the 
microscopic order parameter) is assumed to exist below 
onset, even though long-range order does not. In the 
simple activation model employed here, there is no order 
parameter below onset, since none of the extended states 
in the system are occupied at T = 0. A similar distinc- 
tion appears for the case of the superconductor-insulator 
transition, for which two classes of transitions have been 
observed. For granular films, the order parameter is. non- 
zero at the superconductor-insulator transitionJ23 and 
the superconductor-insulator transition is driven by fluc- 
tuations in the phase of the order parameter. There is 
also evidence for a Bose_glass phase in some homoge- 
neously disordered films ,E3 but other experiments sug- 
gest that the amplitude of the order parameter vanishes 
at the superconductor-insulator transition in this case.Ej 
At the moment, we cannot probe the microscopic order 
parameter in our system, and so we cannot determine 
directly whether or not a condensate exists for coverages 
below the onset of superfluidity. 

C. Thin Superfluid Films 

As discussed above, the rapid change of the chemical 
potential /1 with coverage near n c tends to make quantum 
fluctuations less important in 4 He films than they would 
be if the adsorption energy were less strongly dependent 
on coverage. The large spread in energy scales due to the 
adsorption potential minimizes exchange between atoms 
in the localized layer and those in the overlying superfluid 
film. We believe this lies behind the qualitative success of 
the inert layer model for thicker films, and it is reasonable 
to ask how close to the onset transition one can neglect 
exchange with the non-superfluid film. 

If we ignore the 4 He atoms in the localized layer in 
our aerogel sample, we estimate that the inter-particle 
spacing d for the atoms in the superfluid film is at least 
10 A for T c < 80 mK. This estimate follows from divid- 
ing the number of superfluid atoms by the surface area 
(9.2 m 2 ) of the substrate. Estimates based on the total 
volume occupied by the substrate (0.13 cm 2 ) give an in- 
terparticle spacing on the order of 30 A for a film with 
T c = 80 mKa Both of these estimates indicate that it is 
not unreasonable to look for dilute Bose-gas like behav- 
ior, since the s-wave scattering length characterizing the 
He-He interaction is of order 3 A. 

A full formalism for the crossover from a strongly in- 
teracting superfluid to a dilute Bose gas as the 4 He cov- 
erage is reduced has been developed by Weichman and 
co-workers.E 2 ] The starting point for this theory is the ob- 
servation that thin films of 4 He adsorbed on Vycor show 
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3D XY critical behavior for ..coverages with transition 
temperatures above 120 mK.El Identical critical behav- 
ior (for films with T c 's down 50 mK) has been observed 
in 4 He films adsorbed in 91% porosity aerogel glass.13 In 
both systems, power-law behavior of the superfluid den- 
sity is observed over at least one order of magnitude in 
reduced temperature. Wcichman et al. derive a scal- 
ing function which in principle can be used to test for 
crossover from conventional 3D XY critical behavior to 
dilute Bose gas behavior as the 4 He density decreases. 
The scaling function reduces to the dilute Bose gas form 
Ps(T) ~ p s (0)[l - (T/T c ) 3 / 2 ] as T -> and to the 3D 
XY model form p s (T) ~ p So (l - T/T c ) c (( « 2/3) as 
T — ► T c . Physically, the scaling function accounts for 
the reduction in the size of the 3D XY critical regime 
as T c decreases. Unfortunately, the Bose gas and 3D 
XY limits are sufficiently distinct that a scaling func- 
tion capable of accounting for the crossover can be use> 
to collapse data from essentially any superfluid film 
The two limiting forms, however, suggest another, way 
of testing for crossover to Bose gas-like behavior .Ej The 
quantity lim^^^ p s (T)/(T c — T) diverges in the strongly 
interacting limit but remains finite for the dilute Bose 
gas. In Fig. |(| we plot AP(T)/(T C - T) as a function 
of T/T c for four coverages of 4 He adsorbed on aerogel. 
The dilute Bose gas prediction (assuming the same ra- 
tio AP(0)/T c = 12 psec/mK found for the four films) is 
shown as a dash-dotted line. The divergence near T c is 
pronounced for all of the films. 

There is thus no convincing evidence that 4 He films 
adsorbed in aerogel approach Bose-gas like behavior for 
films with T c 's down to 20 mK. The absence of any ob- 
servable crossover may be due to exchange with the lo- 
calized layer or simply the possibility that our films were 
not sufficiently dilute. Crooker et al. studied films with 
T c 's down to 7 mK in their study of the 4 He- Vycor sys- 
tem and found a progressive decrease in the 3D XY crit- 
ical regime as determined from the same analysis used 
here to produce Fig. |2^. Unfortunately, the superfluid 
signal for a given substrate volume is smaller for aerogel 
than for Vycor, and the thermal time constants are much 
longer in aerogel, so extending our measurements to lower 
coverages and temperatures would not be a trivial task. 
This is unfortunate, since the superfluid transitions in 
4 He-aerogel are "sharper" than those in 4 He- Vycor, in 
the sense that p s goes to zero only 0,2-0.3 mK above T c 
(as determined from power-law fitsu) .as opposed to 1- 
2 mK in Vycor. Kotsubo and Williamso have noted that 
the smearing of the transition in the Vycor case tends 
to mimic crossover to Bose gas-like behavior. There is 
now convincing evidence that the 1-2 mK "foot" in the 
Vycor data was due to macroscopic inhomogeneities in- 
troduced by .-the leaching process used to produce the 
porous glasso A new Vycor study, with more homoge- 
neous samples, would probably be much more conclu- 
sive than the experiment of Crooker et al. in testing 
for crossov.eE to dilute Bose gas behavior. Porous gold 
substrates,H which recent full-pore measurements sug- 



gest are as homogeneous as Vycor, could be an even bet- 
ter substrate for ultra-low temperature measurements. 

The difficulties in interpreting the p s measurements for 
very thin superfluid films have been compounded by the 
absence of a complementary study of other thermody- 
namic quantities in this coverage regime. Our heat ca- 
pacity measurements for He- Vycor begin to address this 
problem, although the lowest transition temperature we 
were able to achieve was only 21 mK, significantly higher 
than the lowest T c (~ 7 mK) reached in the p s study of 
Crooker et al£3 

In Fig [H], we show heat capacity data for two super- 
fluid coverages, n — 28.8 /imoles/m 2 (T c = 85 mK) and 
n = 27.6 /imoles/m 2 (T c = 20.7 mK). The heat capacity 
is divided by temperature and is plotted as a function 
of the reduced temperature (T — T c )/T c . The heat ca- 
pacity for both coverages crosses over from a predomi- 
nantly linear dependence on temperature above T c to a 
quadratic dependence for T < T c . In addition to these 
"backgrounds", there is a peak in the upper data set 
(T c = 85 mK) at the transition. A peak can only barely 
be resolved for the film with T c = 20.7 mK. This-tapd 
is consistent with that observed in earlier studies,! 28 !] 29 ! in 
which the magnitude of the heat capacity peak for 4 He 
films adsorbed in Vycor was shown to scale according to 
the hyperuniversality relation,cfJ 



cm m 3 



c.(t) e(t) 3 



(22) 



for two systems, denoted primed and unprimed, within 
the same universality class. In Eq. ^2], t is the reduced 
temperature (T c — T) /T c and £ is the correlation length. 
Correlation lengths for the 4 He-VycoiE-system can be es- 
timated from the Josephson relation,c3 



kriTm' 



1 



(23) 
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and p s data fr o m the studies of Bishop et aZ.,0 and 
Crooker et aLE3E3 Based on these studies, we estimate 
that the correlation length for a given reduced temper- 
ature will be approximately twice as large for the film 
with T c = 20.7 mK than for the film with T c = 85 mK.El 
According to the hyperuniversality relation Eq. |2^, the 
heat capacity peak for the thinner film should therefore 
be eight times smaller than for the thicker film, which is 
consistent with Fig. |2l]. 

A rigorous test of Eq. |2^ for 4 He films in Vycor would 
require fitting each of the heat capacity data sets to power 
laws near T c . This is impractical at the moment due to 
both the noise in the data and the ambiguities in sub- 
tracting a non-singular background. Although the noise 
could be reduced in a better experiment, the interpreta- 
tion of the background heat capacity poses a more fun- 
damental problem. In Fig. ^l], the data for T c = 20.7 mK 
show a cusp at T c , where C(T) crosses over abruptly from 
quadratic to linear dependence on temperature. This 
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cusp may be the remnant of the fluctuation peak seen 
for thicker films, in which case it should round out as 
the superfluid coverage is further reduced. Another pos- 
sibility is that the heat capacity is singular even in the 
absence of order parameter fluctuations. Although this 
type of behavior would emerge in any mean-field treat- 
ment, for an interacting 4 He film we would expect a jump 
discontinuity in C'(T) in addition to a change in slope. 
(The ideal Bose gas should show only a simple cusp at 
T c , but of a different form than that seen in Fig. gjj .) 

VI. CONCLUSIONS 

The superfluid density and heat capacity measure- 
ments presented in this paper indicate that a 4 He film 
adsorbed on a disordered substrate undergoes a transi- 
tion from a Bose insulator to a superfluid at the critical 
coverage n c . The gap in the insulating phase decreases 
to zero as the onset coverage is approached from below. 
Our data show no evidence of a Bose glass phase. In 
considering thin superfluid films, the rapid change of the 
chemical potential with coverage near n c appears to min- 
imize the role of quantum fluctuations, so that the ordi- 
nary superfluid transition at T c remains the dominant 
critical point. Evidence of critical fluctuations near T c 
is seen even for the thinnest superfluid films studied in 
these experiments. 

Although some of the questions raised in this work 
could be addressed in experiments at even lower tem- 
peratures and smaller reduced coverages, the use of dif- 
ferent substrates may be more rewarding. We have ar- 
gued here that the strong substrate potential in porous 
glasses reduces the importance of quantum fluctuations 
near the critical coverage. Weakly binding alkali metal 
substrates 4 ^ may ameliorate this effect. In principle, the 
use of ordered substrates would allow a test of the clean 
limit of the various dirty boson models, provided that 
the adsorption potential is sufficiently weak. Such a sys- 
tem might be achieved by pre-plating a clean high-surface 
area substrate, such as exfoliated basal-plane graphite, 
with molecular hydrogen. Finally, even on a strong- 
binding ordered substrate such as graphite, the adsorp- 
tion potential changes only weakly within each mono- 
layer, so that quantum effects should be more readily 
observable. The phase-separated superfluids which ex- 
ist within the second and third monolayers of 4 He on 
graphite may therefore be reasonable analogs of granular 
superconductors.cll 
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FIG. 1. Generic phase diagram in the density-temperature 
plane for a 4 He film adsorbed on a disordered substrate. 



FIG. 2. Phase diagram for the Bose-Hubbard model in the 
limit of strong disorder in the plane of chemical potential /i 
and exchange strength J. Both coordinates are scaled by the 
strength V of the repulsive interaction. The transition from 
the Bose-glass (BG) to the superfluid (SF) phase occurs at 
the critical chemical potential /i c as the system travels along 
a path of nearly constant J, as would be appropriate for the 
onset transition in a 4 He film. 



FIG. 3. The resonant period (open circles) and ampli- 
tude (closed circles) of the torsional oscillator as a func- 
tion of temperature during a drift measurement at a cov- 
erage of 35.7 /imoles/m 2 . The period has been offset by 
-3.992359 msec. 



FIG. 4. A cross-section of the 4 He- Vycor calorimeter. All 
pieces are made of silver except as indicated. 



FIG. 5. The CMN magnetization thermometer. 



FIG. 6. The response of the calorimeter following heat 
pulses at 102.1 mK (upper panel) and 18.6 mK (lower panel). 
The solid curves show exponential fits to the decay curve at 
long time scales. Note the deviations at short time scales in 
the 18.6 mK data. 



FIG. 7. Heat capacity, divided by temperature, for a cover- 
age of 27.6 /u,moles/m 2 of 4 He on Vycor. The cusp at 20 mK 
is a signature of the superfluid transition. 



FIG. 8. The superfluid period shift AP as a function of 
temperature for twelve coverages of 4 He adsorbed on 91% 
porosity aerogel glass. 
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FIG. 9. The superfluid transition temperature T c and the 
period shift at T — 0, AP(0), as a function of the the 4 He 
coverage for films adsorbed in aerogel. The solid curves are 
fits to the power laws of Eqs. [L5j and ha 



FIG. 10. The superfluid transition temperature T c as a 
function of the period shift at T = 0, AP(0). The solid curve 
is a fit to Eq. [I]]. 

FIG. 11. AP(0) as a function of the 4 He coverage for films 
with transition temperatures up to 1 K. The solid curve is a 
linear fit to the data for films with T c 's greater than 150 mK. 

FIG. 12. The superfluid transition temperature T c as 
a function of the areal superfluid density at T = 0, 
Ps(0) for films with transition temperatures up to 1 K. 
The solid curve is a linear fit to the data for films 
with T c < 300 mK. The dashed curve is the KTN line 
Ps{T~)/T c = 8.73 pinoles m -2 K _1 . 

FIG. 13. The heat capacity, divided by temperature, is 
shown as a function of temperature for 10 coverages of 4 He 
adsorbed on Vycor glass. The coverages in /imoles/m 2 are 
given in the legend. Non-superfluid coverages are shown in 
the upper panel and superfluid coverages are shown in the 
lower panel. 

FIG. 14. The superfluid transition temperature T c (circles) 
and the crossover temperature Tb (triangles) defined in the 
text are shown as a function of coverage for 4 He on Vycor. 
The onset coverage n c determined by extrapolating the T c 
data is indicated with the arrow. 

FIG. 15. A three-dimensional representation of the C/T 
data for 4 He in Vycor for T < 200 mK over the coverage 
range of our experiment. 



FIG. 18. Low temperature heat capacity data for 
non-superfluid coverages of 4 He in Vycor. The solid curves 
are fits to the activation model discussed in the text. The cov- 
erages in /imoles/m 2 are indicated in the legend. The inset 
shows the gap (see Eq. ^) as a function of the coverage. 

FIG. 19. The coefficients A (circles) and B (triangles) in 
fits of the high-temperature (T > T c or T > 2Tb) heat capac- 
ity data to the form C(T) — AT + BT 2 . The arrow indicates 
the onset coverage n c . 

FIG. 20. Period shift data for four coverages, scaled as de- 
scribed in the text, compared with the predictions for a dilute 
Bose gas. The Bose gas prediction is shown as a dash-dotted 
line for the case AP(0)/T c = 12 psec/mK. The transition 
temperatures for each of the data sets are indicated in the 
legend. 

FIG. 21. The heat capacity divided by temperature shown 
as a function of the reduced temperature (T—T c ) /T c for cover- 
ages of 28.8 /imoles/m 2 (T c = 85.0 mK) and 27.6 /imoles/m 2 
(T c = 20.7 mK). 

TABLE I. Summary of the asymptotic forms for the super- 
fluid density at T=0, p a (0), the superfluid transition temper- 
ature T c , and the heat capacity C(T) in the limits n — > n c 
and T — > 0. The values given are the relevant exponents 
for the Bose glass and inert layer models in two and three 
dimensions. 



Asymptotic Form 


Exponent 


Bose Glass 


Inert Layer 






d=2 


d=3 


d=2 d=3 


p a (0) ~ (n - n c ) c 


c 


> 2 


> 8/3 


1 1 


T c ~ (n - n c ) w 
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> 2 
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1 1 


T c ~[p s (0)f 
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3/4 


1 2/3 
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FIG. 16. The period shift data for coverages of 4 He in aero- 
gel with the period shift scaled by AP(0) and the tempera- 
ture scaled by T c . Open triangles correspond to a coverage 
of 52.7 /imoles/m 2 , for which T c = 955 mK. These data were 
obtained during a separate run. All closed circles are for cov- 
erages with T c < 200 mK. The lowest coverage in this figure 
is 35.2 /imoles/m 2 , for which T c — 26 mK. 



FIG. 17. A model single-particle densiW of states for 
non-superfluid 4 He films adsorbed in Vycor £3 The states be- 
low Eo are localized, and the states above Eo + A are ex- 
tended. 
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TABLE II. Torsional oscillator parameters for the 
4 He-aerogel measurements. Additional details are given in 
Ref. ^. The substrate velocity is estimated from the cell ge- 
ometry and the various amplifier gains. 

Frequency (Hz) 

Moment of Inertia 
(g-cm 2 ) 

Cell Volume (cm 3 ) 

N 2 BET Surface 

Area (m 2 ) 
Mass Sensitivity 

(nsec/pigram) 
Superfluid Mass 
Sensitivity (nsec//igram) 

Q at 10 mK 

Maximum Substrate 
Velocity at 10 mK (cm/sec) 

a Limiting value at high coverages. See Ref. 



TABLE III. Experimental values of the exponents x, and 
w for aerogel and Vycor. Tie Vycor exponents are based on 
fits to the data of CrookerEJ 



Exponent 


Aerogel 


Vycor 


c 


1.64 ±0.04 


2.5 ±0.3 


X 


0.95 ±0.02 


0.85 ±0.02 


w 


1.59 ± 0.06 


1.76 ±0.01 



250 
0.34 
0.13±0.01 
9.2±0.9 

0.211 

0.022 a 
1.6 x 10 6 
2 x 10~ 2 
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